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I. INTRODUCTION 



In a previous paper we presented a formalism that enhances the predictive power of 
perturbative Quantum Chromodynamics (pQCD) with respect to the description of semi- 
inclusive hadroproduction at lepton-hadron colliders. We continue that study here, in 
part because new data has become available from HERA. The reaction that we consider 
is / + A — > I' + B + X, where /, A and I', B are initial- and final-state leptons and hadrons, 
respectively. Calculations in low orders of pQCD are only partially successful in explaining 
the abundant experimental data on semi-inclusive deep-inelastic scattering (sDIS) produced 
in recent years at the DESY-HERA ep collider Satisfactory agreement between the 

theory predictions and the experimental data requires scrupulous and systematic calculation 
of higher-order radiative corrections and, when necessary, all-order summation of poorly con- 
verging perturbative series. A formalism that was originally proposed in || and extended in 
[]J serves the latter purpose. It provides an improved theoretical description of nearly mass- 
less final-state hadrons that escape near the direction of the momentum of the electroweak 
vector boson exchanged between the leptonic and hadronic parts of the sDIS process. In 
other words, this formalism describes production of hadrons in the region of current frag- 
mentation. In this kinematical regime, multiple radiation of soft and collinear partons is a 
dominant factor influencing the behavior of experimental observables. We have shown that 
the effects of multiple parton radiation are most evident in pseudorapidity spectra in the 
center-of-mass frame of the initial hadron and the exchanged electroweak boson. 

The main goal of this paper is to demonstrate that resummation of large logarithmic 
corrections in the current region of sDIS improves the theoretical description of many aspects 
of this process. We compare predictions made within the resummation formalism to various 
sets of experimental data and to predictions obtained for these observables from a next-to- 
leading order (NLO) perturbative QCD calculation. 

In Section[n], we present the equations for the 0(a s ) and resummed cross sections that 
will be compared to experimental data in the subsequent Sections. In Section|ITI|, we compare 
theoretical predictions to new HERA data ||] on the energy flow in the hadronic center-of- 
mass reference frame. In Section |TV|, we discuss the impact of multiple parton radiation on 
various aspects of semi-inclusive charged particle production: in particular, resummation 
effects on transverse momentum distributions and on the correlation between the trans- 
verse and longitudinal momentum components of the observed final-state hadron. Section^ 
contains our conclusions. 

II. SUMMARY OF THE FIXED-ORDER AND RESUMMATION 

CALCULATIONS 

In this paper we use notations and variables that were originally defined in Sections II 
and III of [|J. For the convenience of the reader, we briefly list the relevant definitions in 
Appendix A. 

Following |7|||, we introduce a momentum scale qx related to the space-like vector q!f 
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given by 



t = ?- A ±£5- A JL£i-. (!) 

Pa -Pb Pa- Pb 
The square of the scale q T is defined by 

Qt = -QtQt/x- (2) 

Throughout this paper we will consider observables in the hadronic center-of-mass (hCM) 
frame, in which the four-vector gf has a simple interpretation. This frame is defined by the 
condition p A + q = 0. Let the direction of the z-axis in the hCM frame coincide with the 
direction of q. Then the energy component and z component of q% vanish. The transverse 
components form a two-dimensional vector q^ that is simply the transverse momentum pr 
of the final hadron B rescaled by the Lorentz scalar —z = — (p B ■ Pa)/(q • Pa)', that is, 

gf = {0,g T ,0} = {0,-^,0}. (3) 

z 

The variable z controls the fragmentation of partons in the final state. Thus the vector 
qT is closely related to the transverse momentum pr, except that it is not sensitive to the 
magnitude of the momentum vector ps, i.e., not sensitive to the collinear final-state frag- 
mentation. Rather, qr determines the direction in which the final-state hadron is produced. 
In both the hadron Breit frame and the hCM frame, the magnitude of this vector q T is a 
function of the pseudorapidity of the final-state hadron. In particular, we will later refer to 
the relationship between qr and the pseudorapidity r] cm of the observed hadronic final state 
in the hCM frame, which is 

q T = We^ cm , (4) 

where W 2 = (p A + q) 2 = Q 2 (l/x — 1) is the square of the center-of-mass energy of the 
virtual photon and the proton. The multiple parton radiation that we are discussing here is 
important in the limit f] cm — > +oo, or, equivalently, qT — > 0. 

To see how convergence of the perturbative series breaks down in the small-g^ re- 
gion, we start by examining the hadron-level sDIS cross section doBAl{dxdzdQ 2 dq\-.dip). 
This cross section can be expressed in terms of convolutions of parton-level cross sections 
ddba/ '{dxdzdQ 1 [ dq T d(p), and phenomenological long-distance hadronic functions: 

do-BA _ ^ f 1 d£ b ^ f 1 d£g p , > da ba (^ F ) 

dxdzdQ 2 dq 2 T d V ~^Jz & B/b[ ^ M L i a alAKa ' dxdzdQ 2 dq 2 d V ' 1 J 

Here F a / A (£, a , A*f) denotes the parton distribution function (PDF) for a parton a in the hadron 
A, and Ds/biCb, A*f) is the fragmentation function (FF) for the hadron B that is produced in 
the fragmentation of a parton b. The notation /xp is used for the scale of factorization. The 
sum over the labels a,b includes contributions from all parton types, i.e., g,u,u,d,d, . . . . 
In the following, a sum over the indices i,j will include contributions from active flavors of 
quarks and antiquarks only, i.e., it will not include a gluonic contribution. 
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Both the hadron-level and parton-level cross sections can be expanded into sums over 
leptonic angular functions A p (i/),ip) (which are listed in Appendix A), as 



(Loba 



1 



dxdzdQ 2 dq T d(p 

dd ha 
dxdzdQ 2 dq T dip 



Y^ P VBA(x,z,Q 2 ,q 2 T )A p (^,v), (6) 



1 



£ p Vu>(x,z,Cf,<£)A p (i/;,<p). (7) 
P =i 



The coefficients p Vba (or p Vb a ) of the angular functions A p (ip, <p) are independent of one 
another. If <C 1, but qx ^ 0, the perturbative parton-level cross section is dominated 

by the term with p — 1. The structure function for this term l Vba(x, 2, Q 2 , q T ) behaves as 
l/q T times a series in powers of a s and logarithms ln(gfyQ 2 ), 

The constants o"o and i 7 }, which were introduced in JT], are 

(*)■ < 9 > 

(10) 



O"0 



A7rS eA x 2 V 2 
2 Q 2 



The coefficients (x, z) are generalized functions of the variables x and 2. 

Convergence of the series in (||) deteriorates rapidly as qx/Q — > 0, because of the growth 
of the terms (q^ 2 ) In" 1 (q^ / Q 2 ) . As a result, the structure function 1 Vb a (x, z, Q 2 , q?) has a 
non-integrable singularity at qr = 0. Its asymptotic behavior is very different from that of 
the structure functions 2 ' 3,4 H a (x, z, Q 2 , gf*), which are less singular and, in fact, integrable at 
qx = 0. This singular behavior of 1 Vfe a (x, z, Q 2 , q?) is generated by infrared singularities of 
the perturbative cross section that are located at qx = 0. These singularities must be either 
explicitly canceled or factored into the long-distance hadronic functions in order to obtain a 
finite cross section. 

There exist two qualitatively different approaches for handling such singularities. The 
first approach deals with the singularities order by order in perturbation theory; the second 
approach identifies and sums the most singular terms in all orders of the perturbative expan- 
sion. We will begin by discussing the first approach, in which singularities are regularized 
independently at each order of the series in a s . In the crudest approximation, one keeps only 
terms of the order 0(a s ) and drops all higher-order terms. The singularity in the 0(a s ) 
part of the asymptotic expansion @ can be regularized by introducing a "separation scale" 
q^ p and considering the fixed-order cross section separately in the regions < qx < q'x P an d 
It > 1x P ■ The value of q^ v should be small enough for the approximation (^) to be valid 
over the whole range qx < qx V ■ 

In the region < qx < qx % \ we can apply the MS factorization scheme to take care of the 
singularities at qx = 0. In the MS scheme, the regularization is done through continuation 
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of the parton-level cross section to n = 4 — 2s dimensions. The n-dimensional expression for 
the asymptotic expansion (||) of 1 Vt, a (2, z, Q 2 , q^) is 



27^,4 V z 

The scale parameter \x in ([11]) is introduced to restore the correct dimensionality of the 
parton-level cross section da ba /(dxdzdQ 2 dq^ T ~ 2 dip). The soft and collinear singularities ap- 
pear as terms proportional to l/e n when n — > 4. The soft singularity in the real emis- 
sion corrections cancels with the soft singularity in the virtual corrections. The remaining 
collinear singularities are absorbed into the partonic PDFs and FFs, which therefore must 
be redefined at each order of the perturbative calculation. When the partonic PDFs and 
FFs are subtracted from the partonic cross section da, the remainder is finite and inde- 
pendent of the types of the external hadrons. We denote this finite remainder as (da)hard- 
The convolution of (da) hard with the hadronic PDFs and FFs yields a cross section for the 
external hadronic states A and B. There is a certain freedom in the precise definition of 
the collinear contributions, so that the "hard" part depends on an arbitrary factorization 
scale fip through terms like P ab ln(fi F /K), where P a b(0 are 0(a s ) splitting functions ||, 
and K is some momentum scale in the process.^] The dependence on the factorization scale 
in the hard part is compensated, up to higher-order terms in a s , by scale dependence of the 
long-distance hadronic functions.^ 

After the cancellation of soft singularities and factorization of collinear singularities, one 
can calculate analytically the integral of ^Vbajhard over the region < q F < q^ P - At 0(a s ) 
this integral is given by 

r(q s T ep ) 2 , ^ x a Fi „ 2 f x ^ LO aJM uJnlo\ n 9 n 

J ^ { V H hard - ^7 Yi { V ba,J + V ba,j ^ (12) 

where j — u,u,d,d, . . . ; the LO and NLO structure functions are 

l V h L a °i = 5(l-z)5(l-x)8 bj 5 ja , (13) 



lf> NLO 
V ba,j 



C F In 2 j^ijy - 3C F In j^py J K 1 ~ ~ x)S bj 5 ja 

,2 



+ In — 6(1 - z)5 bj P ja (x) + P hj (z)8(l - x)8 ja 
[zq T J \ 

+ 5(1 - z)5 bj cj} l) (x) + c°f (1) (z)5(l - x)5 ja . (14) 



1 The scales [ip and [i are related as /sp = Aire 7B /i 2 , where ~/e is the Euler constant. 

2 In our calculations we do not distinguish between factorization scales for PDFs and FFs nor 
between the factorization and renormalization scales. 



5 



The coefficient functions c^' out( - 1 ^(^) that appear in x V h N a L ^° are given by 



1(1-0-2(5(1-0 



(15) 



4 )m (0 = ^(i-0, (16) 

4f"(0 = (17) 

Now consider the kinematical region q<r > <?r P , where the approximation ( |TTD no longer 
holds. In this region, ( 1 Vb a )/ iar d should be obtained from the exact NLO result. With this 
prescription, the integral over q\ can be calculated as 



1 



max " 2 



It 2 da ba 
aqf 



dxdzdQ 2 dq^dcp 

(9t P ) 2 o /t- \ /-maxcj 2 



r\q T ) / ^ \ /"maxqj, / ^ \ 

A 1 (i>,i P ){ dql( l V ba ) h + / dq 2 T ( l V ba ). 

'{JO V 'hard J(q£ p )a ^ V Jl 

j finaxgi, . ^ . 

+ E^^) / ^( P Ha) 




p=2 

where maxg^ is the maximal value of q\ allowed by kinematics. The first integral on the 
right-hand side is calculated analytically, using the approximation (0); the second and 
third integrals are calculated numerically, using the complete perturbative result of the 
order 0(a s ). The numerical calculation is done with the help of a Monte Carlo integration 
program written in the style of the package ResBos used earlier for resummation in vector 
boson production at hadron-hadron colliders [|T7J . 



A significant failure of the calculational procedure in ( jig) is that it does not reproduce 
the correct shape of g^-dependent differential cross sections in the region q? <C Q. The 
fundamental problem is that the terms in (|11|) with small powers of a s do not reliably 
approximate the complete sum in the region qr ^ Q, where the large magnitudes of the 
logarithmic terms \n m {q^/Q 2 ) quickly invalidate the expansion in powers of a s . In other 
words, a truncated sum in powers of a s (cf. Eq. (|T2|) ) cannot be a reliable prediction for the 
observable cross section in the region qr/Q "C 1. 

The problem mentioned in the previous paragraph justifies the second approach to the 
regularization of the singularities at qx = 0, in which large logarithms in (|ll|) are summed to 
all orders. According to the formalism developed in |H|]I^fl and adapted to the sDIS process 
in m, a better approximation for this sum is provided by the Fourier transform of a 6-space 
function W(b, Q, x, z, Ci, C2, Hf), which sums the dominant terms in (j|) through all orders 
of a s \ 

1 Ha(x,2,g 2 ,g^) = -^-e^W^Q^C^^). (19) 

Here b is a vector conjugate to qr, and b denotes the magnitude of b. 
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From first principles jll|,|l2| , the general structure of Wb a (b, Q, x, z, Ci, C 2 , Hf) is known 
to have the form 



W ba (b,Q,x,z,Ci,C 2 ,iiF) = 

e 2 3 C° b f(z, &, C 2 , K, fi F )e- sP ^ Q ' c ^e- sNP C; n a (x, C u C 2 , K, fx F ), (20) 

3 

where j = u, u, d, d, . . . . The meanings of various factors in Eq. (|20| ) were discussed in 
detail in M. Here we should emphasize that at short distances, i.e., in the limit 5 — * 0, 
the function Wba(b,Q,x,z,Ci,C 2 ,fXF) is dominated by perturbatively calculable functions 
Cja(x, C\, C 2 , b, ftp), ^^(z, Ci, C 2 , b, fxp), and the perturbative part of the Sudakov factor 
S p (b,Q,C\,C 2 ). The 0(a s ) expressions for these functions are listed in Appendix B. 

In the region b ^ b max , where b max can be chosen to be 0.5 GeV -1 , nonperturba- 
tive contributions to Wf, a become important. These nonperturbative contributions are ab- 
sorbed into a function S (b,Q, . . .), which vanishes in the limit b — > but dominates 
Wba(b, Q, x, z, Ci, C 2 , ix p) for b > b max . Until a nonperturbative theory of strong interactions 
is developed, the explicit parametrization of S NP (b, Q, . . .) can only be found from compar- 
ison with experimental data. One possible realization of a smooth transition in the region 
b ~ b max between the perturbative and nonperturbative parts of Wb a (b, Q, x, z, C\, C 2 , Lip) 
is to replace the parameter b in the C-functions and the perturbative Sudakov factor S p by 
the variable 

K = b (21) 

yjl + {b/b max ) 2 

as written in ( p0|) . 

There is a subtle difference between C m and C out . While the initial-state coefficient func- 
tions (x, Ci, C 2 , b, n f) given in Eqs. (|70|) and (ffTD depend on the factorization scale llf 
through a factor In [//p&/& ], the final-state functions {z, C\, C 2 , b, lxf) given in Eqs. ( [72]) 

and ( [75| ) depend instead on In [yLpb/ {b^z)]. The additional term oc ln£ in the functions 
Cfa (z, Ci, C 2 , b, up) becomes large and negative when z — > 0, so that it can significantly 
influence the 0(a s ) contribution at small values of z. As a result, the resummed total rate 
tends to be lower than its fixed-order counterpart for z ^ 0.1. We will discuss this issue in 
more detail in Section|TV|. Similarly, the 0(a s ) part of the NLO structure function l V b ^° 
in fllij) depends on Lip through a logarithm In [n F /(zq^ p ) 2 }. 

The appearance of the additional terms oc In z in the functions C° b ; t{1) and l V b N a f reflects 
the specifics of separation of the 0(a s ) "hard" cross section (a) hard from the collinear contri- 
butions to the FFs in the MS factorization scheme. The easiest way to see the specific origin 
of the In z terms is to notice that the dependence on the parameter fx in the n-dimensional 
expression (|TTD for 1 Vfea(^, z, Q 2 , q^) comes through a factor {2'Kfx/z) A ~ n , rather than through 
a more conventional (2nfx) 4l ~ n . It is this factor that generates the /z^-dependent logarithmic 
terms In \fxpbj (b z)] in the functions Cl^ t<yl \z, C\, C 2 , b,/ip) and l V b ^° . The Cj" -functions 
do not include In z because they are evaluated along the direction z — 1 in the phase space. In 
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contrast, nothing forbids such a term in the functions C b j , in which z can be anything be- 
tween z and 1. This contribution produces terms of the form Pj,j In z in (z, Ci, C2, b, 
which suppress the resummed rate for the production of hadrons with very small energies. 

In Ref. , the terms proportional to In z were left out of the expressions for the C^- 
functions. This error affected the value of the C ou *-function for the energy flow in Ref. |]|; 
the correct expression is Eq. ( j30f) in this paper. We have found that theoretical distributions 
that were shown in the Figures of Ref. [Ij] will be approximately reproduced with the correct 
C£a* ^ the following nonperturbative Sudakov factor is used: 

S? P (6) = In ^ + 6 2 (-4.58 + . (22) 

The explicit parametrization for the coefficient g^^b) * s gi ven m 0- I 11 the next Section, we 
present an updated analysis of the transverse energy flow that is compared to the newer data 
from the HI Collaboration. We will present a new parametrization of the nonperturbative 
Sudakov factor which is in better agreement with the whole set of experimental data. 



III. RESUMMATION FOR THE TRANSVERSE ENERGY FLOW 

As was discussed at length in |]T[, knowledge of the resummed semi-inclusive DIS cross 
section leads immediately to predictions for the pseudorapidity spectrum of the transverse 
energy flow in the hadron Breit frame or the hCM frame. It is advantageous to study the 
energy flows, because they are less dependent on the specifics of final-state fragmentation of 
the scattered partons into the observed hadrons. We therefore start the presentation of our 
numerical results with the comparison of the resummation formalism to the experimentally 
measured pseudorapidity distributions for the transverse energy flow in the hCM frame. 

According to Eq. , the scale for massless particles is simply related to the pseudora- 
pidity f] cm in the hCM frame, if the energy of the photon-proton collision is fixed. Therefore 
the pseudorapidity distribution for the transverse energy flow in the hCM frame d{Ej) / ' dr\ cm 
can be converted into a g-r-distribution d{ET)/dqT by a simple change of variables. From 
the viewpoint of theoretical analysis, it is preferable to convert d(ET)/dr] crn into the qr 
distribution of the "z-flow" [Ol that is defined as 



-22 z dz - ( 23 ) 



dx dQ 2 dqr # Jz min dx dz dQ 2 dqr 

Under the approximation that all particles are massless, the distribution d(Ex) jdif m is 
simply related to the distribution dT, z /dqx, by 

= 4 (24 ) 
drjcm dq T 

In [0 the predictions of the resummation formalism were compared to data on transverse 
energy flow (Et) published in 0. In a recent publication ||, the HI Collaboration at the 
DESY-HERA ep collider presented new data on the distributions d{Ej) /df] cm in the same 
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bins of x and Q 2 as in || and, in addition, in several bins at significantly higher x and Q 2 . 
In this paper, we repeat the analysis of to produce distributions dH z /dqT for the bins of 
x and Q 2 discussed in Q. The procedure of our numerical analysis is identical to the one 
used in [p]]. Our calculations use the parameters of the HERA ep collider. The energies of 
the proton and electron beams are taken to be 820 GeV and 27.5 GeV, respectively. The 
current numerical results are obtained using CTEQ5M1 parton distribution functions . 

By the factorization theorems of QCD, the hadron-level z-flow S 2 can be written as the 
convolution of a parton-level z-flow H z with the PDFs, 



(l.nl()-<l(j-j,d(p ^ Ix £ a a//j4 ^ a ' dxdQ 2 dq\dtp ^ ^ 

Similarly to the sDIS cross section, the z-flow can be expanded in a sum over the leptonic 
angular functions A p (ijj,<p): 



p=i 



(26) 



where the structure functions p V za (x,Q 2 ,q^) for the z-flow are related to the structure 
functions p Vf, a (x, z, Q 2 , q?) for the sDIS cross section by 

p V za (x, Q\ q 2 T ) = £ T zdz PVtaix, z, Q 2 , q 2 ). (27) 
b Jo 

The resummed z-flow is calculated as 

H 

dxdQ 2 dq^d(p S eA 2 J (27r) 2 
where 

W z (b, x,Q)=J2 e 2 f° ut e- s ^Q>*) (C% o F a/A ) (x, K, (jl f ). (29) 



/ 7^^ W M^Q) + ^, (2f 



The functions C 1 -^ in (|29|) are the same as in (^0|). The coefficient C° ut is 



(~ = 1+*C,( * £_ h »£^, (30) 
7T v 16 3 C 2 6o ' 

where b = 2e~ lE and 7^ = 0.577. . . is the Euler constant. The parameter given by 
(pl|) with b max = 0.5 GeV -1 , is introduced in (^) to smoothly turn off the perturbative 
dynamics when b exceeds b max . The numerical results for the resummed z-flow are obtained 
using Ci/C 2 = l^Fb = b . 

The term Y z in (|28| ) is the difference between the complete fixed-order expression at 
0(a s ) for dY Jz j{dxdQ 2 dq\di^) and its most singular part in the limit qx — > 0; that is, 

y^^J^-(^^r-) ■ (3D 



dxdQ 2 dq\dtp \dxdQ 2 dq^dip 



sing 
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The singular part calculated to 0(a s ) is 

\dxdQ 2 dq^d(f j sing S eA vr 2q\ 2tt 



x 



{(Pqq O F j/A )(x, Hp) + (P qg O F p/A ) (x, 



+ 2F iM (a;, ^){c F In % - ^-C F \ + 0((—) \ q 2 T 



(32) 



The formula (|2q) is expected to provide a smooth matching between the resummed z-flow in 
the region <C Q an d the fixed-order 2-flow in the region g^ > Q. However, we find that the 
0(a s ) fixed-order z-flow significantly undershoots the data in the region qr/Q ~ 1. The re- 
summed z-flow shows better agreement with the data, over a wide range of qr/Q, but also lies 
below the data if qr/Q significantly exceeds unity. The discrepancy between the theoretical 
z-flows and the data at qx/Q ~ 1 is illustrated for the bin (Q 2 ) = 33.2 GeV 2 , (x) = 0.0047 
in Fig. and for the bin (Q 2 ) = 617 GeV 2 , (x) = 0.026 in Fig.|. The data at sufficiently 
high qx lies above the 0(a s ) z-flow and the resummed z-flow in both bins. This difference 
signals that calculation of higher order radiative corrections is required to achieve agreement 
between the theory and the data for qr/Q ~ 1- The better consistency between the re- 
summed z-flow and the data suggests that we should use the resummed z-flow up to values 
of qr/Q ~ 1, ie., without switching to the fixed-order expression. We have followed this 
procedure in the derivation of our numerical results. 

We have also found that the agreement between the resummed z-flow and the data in 
the region qr/Q ~ 1 can be improved by making an additional correction that accounts 
for the reduction of the allowed phase space with the increase of qx- This correction was 
discussed in detail in |TJ. Its importance can be understood by noting the steep rise of the 
PDFs as x decreases, at values of x relevant to the HERA measurements. This steep rise in 
the small x region, together with a mismatch between the phase space in the 6-space integral 



of the resummed z-flow (p8 ) and that in the exact fixed-order formula, may produce large 
numerical effects if q? is of the order Q or larger. To correct for this mismatch between the 
resummation and the perturbative phase spaces, x = Q 2 / (2p A ■ q) is replaced in Eqs. ( p9[ ) 
and (132]) by the g^-dependent quantity 

x = — x. (33) 

The replacement x — > x preserves the leading-order kinematics when qx approaches zero; at 
the same time, it simulates the reduction of the allowed phase space when q? becomes com- 
parable with Q. It affects only the convolution integrals (Pb a °Fa/A) an d (Cb a °F a /A), in which 
the gr-independent lower limit of integration x is replaced by the variable x. The result of 
this replacement is that the convolution integrals at large q F do not include large contri- 
butions from very small values of x, which helps to improve the accuracy of the matching 
procedure at qr ~ Q. 

Before turning to the discussion of the experimental data, let us comment on the structure 
of the Sudakov factor S z (b, Q, x) in (|2"P|). As in the case of the sDIS cross section, the Sudakov 
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factor for the z-fiow consists of both perturbative and nonperturbative parts, 

/■CfQ 2 ^772 / n 2 Q 2 \ 

S g (b,Q,x) = I 1 A(a s (-p), Ci) In + B(a s (jl), C u C 2 ) 

+ Sf p (6,Q,x). (34) 

A realistic parametrization for the nonperturbative Sudakov factor S^ p (b,Q,x) can be ob- 
tained by comparison with experimental data, especially with the measured pseudorapidity 
distributions at Q rs 3 — 20 GeV. As Q increases, we expect the data to become less sen- 
sitive to nonperturbative effects and to be dominated by perturbatively calculable parton 
radiation. In [|TJ we found a parametrization of S^ p (b,Q,x) that was consistent with the 
data published in ||. But, because of the limited scope of those data, that parametrization 
was expected to be valid only for 10~ 4 < x < 1CT 2 . Also, the Q dependence of S^ p (b, Q, x) 
could not be tested, since all the pseudorapidity distributions in Ref. || were presented in 
a small range of Q ~ 2 — 6 GeV. Hence, the parametrization of S^ p (b,Q,x) in [p]] used a 
model prediction for the coefficient of In Q/Qo based on the assumption of crossing symmetry 
between the nonperturbative Sudakov factors in sDIS, e + e~ hadroproduction, and Drell-Yan 
processes. 

In its more recent publication Q, the HI Collaboration presented pseudorapidity distri- 
butions of the transverse energy flow for Q 2 up to 2200 GeV 2 . However, as just mentioned, 
the data at such high Q 2 should be rather insensitive to the nonperturbative dynamics, 
e.g., because the nonperturbative effects would be important only at very large values of 
f] cm , for which the resolution of the detector is poor. Thus the data at very high Q 2 is 
not very informative about the Q 2 -dependence of S^ p (b,Q,x) either. To learn about the 
dependence on Q 2 , one needs to look at data in the range Q ~ 5 — 20 GeV. Fortunately, 
the HI Collaboration presented distributions in two bins at intermediate values of (Q 2 ), 
namely (Q 2 ) = 59.4 GeV 2 and (Q 2 ) = 70.2 GeV 2 . Together with the data from Refs. |],f§ 
at lower values of Q, these distributions provide the first direct tests of the Q 2 -dependence 
of S^ p (b,Q,x). We therefore present a numerical value for the coefficient of liaQ/Q that 
yields reasonable agreement with all of the analyzed data. The resulting value of this coef- 
ficient differs significantly from the model parametrization in |]TJ. However, we should not 
draw too strong a conclusion from this difference, because it might be caused by ambiguities 
in the separation of Q 2 dependence and x dependence in the existing data. To draw a strong 
conclusion about the crossing symmetry model, experimental pseudorapidity distributions 
in a larger range of x at intermediate values of Q 2 will be needed. 

From our analysis, we find that the data from Refs. ||f|] are consistent with the following 
parametrization of the nonperturbative Sudakov factor: 

S? P (b, Q, x) = b 2 ( O.O ^ 1 "^ 3 + 0.19 In ^- + cV (35) 
\ x Q J 

where the parameter Q is fixed to be 2 GeV, and where we set C = for reasons explained 
later. In Figs.[l| and |2] we compare the resummed z-flow, calculated according to (^) and 
using the parametrization in (|35]), to the data from M. The z-flow is normalized by the 
inclusive DIS cross section da j (dxdQ 2 ) (which will be denoted by the brief notation <r tot ). 
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Following the recent HI publication Q , we organize these data into two sets of bins for fixed 
x and Q 2 . One set covers 3.2 < (Q 2 ) < 70.2 GeV 2 and 8 x 1(T 5 < (x) < 7 x 10~ 3 , out of which 
we consider the bins with (Q 2 ) > 13.1 GeV 2 . The other set covers 175 < (Q 2 ) < 2200 GeV 2 
and 0.0043 < (x) < 0.11. The transformation of the experimental d(Ej-) / drj cm distributions 
into dH z /dqT distributions and the calculation of the theoretical distributions dYi z /dqT are 
done using the mean values (x) and (Q 2 ) for each bin published in |4j]. In seven of the nine 
low-Q 2 bins (Fig. [I]) we also show the earlier published data from [0]. 

Let us discuss the features of the data presented in Figs.|l| and ||. First, the data in the 
low-Q 2 bins is significantly influenced by nonperturbative effects and therefore is sensitive 
to the details of the parametrization of the nonperturbative Sudakov factor S^ p (b,Q,x). 
This feature can be seen from the abundance of data points around the maximum of the qx- 
distribution, where the shape is mainly determined by the nonperturbative Sudakov factor. 
Also, the low-Q 2 data from HERA is characterized by small values of x, between 10~ 4 and 
10 -2 . As we noticed in our earlier work |T|], for the theory to be consistent with the data in 
this range of x, the nonperturbative Sudakov factor must increase rapidly as x — > 0, at least 
as l/y/x (cf. Eq. (|22|)). In our new analysis, we find that growth of S^ p (b,Q,x) as 1/x at 
small x is in better agreement with the more recent data from Q. 

Secondly, the data in the high-Q 2 bins of Fig. Q shows a behavior that is qualitatively 
different from Fig. |TJ. In the region covered by the experimental data points, the distri- 
bution is a monotonically decreasing function of qT, which shows good agreement with the 
resummed z-flow over a significant rangeQ of qT- In the region < 10 GeV, i.e., where 
the maximum of the qx distribution is located and where nonperturbative effects are impor- 
tant, the experimental gx-bins are too large to provide any information about the shape of 
dYizjdqrp. Thus, as mentioned earlier, the published high-Q 2 z-flow data from Ref. 0] is not 
sensitive to the dynamics described by the nonperturbative Sudakov factor S^ p (b,Q,x). 

A third comment is that most of the high-Q 2 data points in Fig. ^ correspond to (x) > 
10~ 2 . If the resolution of the HI measurements at high Q 2 were better in the small-g^ 
region, then the high-Q 2 data would also reveal the behavior of S^ p (b, Q,x) at large x. But, 
as mentioned before, the published data in the high-Q 2 bins are not very sensitive to the 
shape of the z-flow at small g^. Therefore it is not possible to impose any constraints on 
S^ p (b,Q,x) at large values of x, except that it should be positive, S^ p (b,Q,x) > 0. For 
this reason we have chosen a parametrization of the x-dependent part of S^ p (b, Q,x) in (|3~5]) 
such that S^ p (b,Q,x) grows approximately as 1/x as x — > and is positive for all x. For 
the same reason, we chose C = in the parametrization fl3"5|). Although the most general 
parametrization of S^ p (b, Q, x) can have C^O, the current data cannot distinguish between 
the parametrization ( |55D with C = and C ^ 0, as long as the value of C is not very large. 



3 We point out once again that both the 0(a s ) and resummed z-flow lie below the data at very 
large qx, in all bins of x and Q 2 in Figures [l] and 0. 
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How trustworthy is the resummed z-flow at large q^l 



As noted earlier, the 0(a s ) fixed-order z-flow is much larger than the data in the region 
Qt/Q "C 1 and much smaller than the data in the region q T /Q > 1. In the small-g^ region, 
the resummed z- flow is, by its construction, more reliable than the fixed-order result. In the 
large-gy region, the resummed z-flow, with the kinematical correction (|33|) included, is also 
in better agreement with the data than the fixed-order calculation. But theoretically, the 
resummed z-flow at large qr/Q is not absolutely trustworthy, because it does not include 
those parts of the fixed-order z-flow that are sub leading in the limit q? — > 0, but which might 
be important at large q?- If the NLO result were in a good agreement with the data at large 
qx, it would be justified to consider it a more reliable prediction in this region. But since 
the 0(a s ) contribution is systematically smaller than the data, higher-order corrections are 
presumably necessary in order to describe the region ~ Q reliably. 

A systematic approach for improving the theoretical description of the large-g-r region 
would require inclusion of the complete 0{a 2 s ) terms in both the fixed-order and resummed 
z-flows. But because such a calculation is not available, it might be beneficial to use the 
resummed z-flow as a better theoretical prediction both in the region qr/Q "C 1, where 
application of the resummation formalism is fully justified, and for q^ up to several units of 
Q, where the resummed z-flow agrees with the data better than the fixed-order one. Then 
the use of the resummed (^-distributions of the z-flow will provide more reliable predictions 
for other observables relevant to the sDIS process. 

As an example, resummation can improve the reliability of the theoretical prediction 
for the azimuthal asymmetry of the z-flow. The 6-space resummation formalism affects 
only the coefficient 1 V z a of the angular function Ai(ij),ip). This coefficient is the one that 
dominates the yj-integrated z-flow in the small-g-r region, where the energy flow is the most 
intense. On the other hand, the main goal of the measurement of angular asymmetries 
is to study structure functions other than 1 V z a, e.g., those corresponding to the angular 
functions As(tp,(p) = — cos (p sinh 2ip and A^ip^ip) = cos2<y9sinh 2 ?/;. By using a better 
approximation for the coefficient 1 V z a, it is possible to measure the coefficients 3 ' a V z a more 
reliably. Conversely, by knowing that the all-order resummation effects are important in 
the region of small qx and by concentrating on the region where qx is of the order Q or 
larger, one may find angular asymmetries that are well approximated in the lowest orders 
of perturbative QCD. We discuss the impact of resummation on the angular asymmetries in 
more detail in another paper \T%\. 



IV. DISTRIBUTIONS OF CHARGED PARTICLE MULTIPLICITY 

Let us now turn to the discussion of semi-inclusive DIS cross sections. Although the re- 
summation formalism, as outlined in Section |J, can describe the cross section for any massless 
final-state particle (provided that the fragmentation functions for this particle are known), in 
this Section we will concentrate on distributions of the charged particle multiplicity, defined 

as 

1 da (A + I — h± + V + X) 
^ d& ' (36) 
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Here G is some kinematical variable, such as the variable q T in Eq. (|2|), the transverse mo- 
mentum px of the final-state charged particle in the hCM frame, or the Feynman variable 
xp, defined by 

In (|37|) p z B is the longitudinal component of the momentum of the final-state hadron in the 
hCM frame. For small values of q F , i.e., in the region with the highest rate, 

x F ss z. (38) 

In our calculation we assume that the charged particles registered in the detector are 
mostly charged pions, kaons and protons. Therefore the cross section for charged particle 
production can be calculated using (|5]) with the replacement of the fragmentation functions 
D B /b(€b, Vf) by 

D h ±/b(£b, Hf) = D B/b(£b, A*f)- (39) 

B 

The fragmentation functions D B / b (^j„ fj,p) are known reasonably well only for > 0.05 — 0.1 
|l9|-pT|. Thus, the formalism presented here is applicable to the production of charged 
particles with sufficiently large energies, i.e., for z > 0.05. 

Certain experimental distributions are readily available from the literature |5jl7,18|, such 



as da/dpT, da/dxF, as well as distributions for the average transverse momentum (p T ). 
However, the "experimental" qp distributions must currently be derived from pseudorapid- 
ity distributions by using Eq. (|). Although the distributions da/dpp and (p T ) are quite 
sensitive to resummation effects, they cannot be interpreted as easily as the distributions 
da/dqp, primarily because the distributions da/dpx and (jp T ) mix resummation effects at 
small values of qp with perturbative contributions from the region qp/Q ~ 1. The most 
straightforward way to study the effects of multiple parton radiation would be to consider 
the qp (or pseudorapidity) distributions that satisfy the additional requirement z > 0.05 — 0.1 
and that are organized into small bins of Q 2 and x. Unfortunately, such distributions have 
not been published yet. 

Because the experimental qp distributions are unavailable, we have decided to undertake 
a simpler analysis than the one presented for the energy flow. Our goal here is to understand 
how the multiple parton radiation could affect various aspects of charged particle production. 
For this purpose we focus our attention on data from the ZEUS Collaboration ||, which 
presents the charged particle multiplicity in a phase-space region characterized by the mean 
values (W) = 120 GeV, (Q 2 ) = 28 GeV 2 , and the additional constraint z > 0.05. These 
values of (W) and (Q 2 ) translate into an average value of x = 1.94 x 10~ 3 . Within the 
framework of a simple model for nonperturbative effects at small qx, we will demonstrate 
that resummation describes qualitative features of this set of experimental data better than 
the fixed-order calculation. In all of the cases presented, our strategy is to compare the 
resummed multiplicity to that from the next-to-leading order calculation. In the numerical 
analysis, the multiplicity was calculated using the CTEQ5M1 PDFs [|l6j and the FFs from 
For the resummed multiplicity, the "canonical" combination G\jCi = fiFO = 2e _TB 
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was used. Unless specified, the NLO cross section was calculated according to (0), with the 
parameter q^ v = 1.2 GeV, and for the factorization scale \ip = Q. 

As in the case of the z-flow, the resummed charged particle multiplicity may suffer from 
matching ambiguities at qr/Q ~ 1- In Section|ITI|, we found that the resummed z-flow is 
in better agreement with the experimental da/dqT distributions than the NLO z-flow, for 
the whole range qr/Q ~ 2 — 4. That result suggests that it might be preferable to use the 
resummed z-flow in the whole range g-r/Q ~ 2 — 4 as a better theoretical prediction, until 
the 0{o?g) prediction for the z-flow in the region qp/Q ~ 1 becomes available. In the case of 
the charged particle multiplicity, the resummed cross section, which is calculated according 
to the formula 

dxdzdQ 2 dq 2 d^ = -S^^^J W? 6 WsA{b, x, z, Q) + Yba, (40) 

overestimates the experimentally measured rate for the production of charged particles with 
Pt > 2 GeV. This discrepancy indicates that the resummed cross section in the region 
Qt/Q ~ 1 is too high, so that switching to the perturbative cross section in this region is in 
fact required. Therefore, we have chosen to use the resummed cross section for qp < 5 GeV 
and switch to the next-to-leading order cross section for qp > 5 GeV. 

As in the case of the z-flow, the shape of the qp distribution for the charged parti- 
cle multiplicity at small values of qp depends strongly on the unknown nonperturbative 
Sudakov factor S NP (b, Q, x, z). For the purposes of this study, we introduce a preliminary 
parametrization of the nonperturbative Sudakov factor for the fixed values of x = 1.94 x 10~ 3 
and Q 2 = 28 GeV 2 , i.e., the values that coincide with the average values of x and Q 2 in ||. 
This z-dependent parametrization is 

S NP (b, Q 2 = 28 GeV 2 , x = 1.94 x 1(T 3 , z) = b 2 ^0.18 + 0.8 (1 ~f^ . (41) 

Since the ZEUS Collaboration did not publish pseudorapidity distributions for the charged 
particle multiplicity i\j o to1 )da j dif m in bins of varying z, we had to deduce information 
about the z-dependence of S NP from the less direct distribution of (p 2 ^) vs. xp presented in 
Fig. 3c of ||. This distribution, known gull" for its characteristic shape (Fig. |3]a), 

can be converted into the more illustrative distribution of (q 2 ^) vs. xp (Fig. §b). Since the 
major portion of the registered events comes from the region q\jW 2 1, or xp w z, a first 
estimate of the experimental data points for the distribution of (q^) vs. xp can be obtained 
by assuming that 



1t) (z) 2 (x f )*' 



where (xp) denotes central values of xp in each bin in Fig. |3]a|] We refer to the resulting 
values as "derived data". 



In principle, a more accurate experimental distribution (qj,) vs. xp can be determined by its 
direct measurement. 
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Note that the shapes of (pp) vs. xf and (qp) vs. xf are quite different. The transforma- 
tion from Fig.|^a to Fig.|3|b shows immediately that the wing-like shape of the distribution 
of (pp) vs. xp should be attributed to a purely kinematical effect, namely an extra factor 
1/z 2 which is absent in the distribution of (q^) vs. xp. Once this extra factor is removed, 
we see from Fig. ^b that (qp) increases monotonically and rapidly as z approaches zero. In 
other words, the qp distribution broadens rapidly when z decreases. This behavior is approx- 
imately realized by the simple z-dependent nonperturbative Sudakov factor S NP (b,Q,x, z) 
given in Eq. (f4T|). 

The shape of the resummed cross section at qp < 5 GeV and the parametrization of 
S NP (b,Q,x, z) were chosen to maximize the agreement between the resummed distribution 
of (qp) vs. xf and the "derived data" (c/. Fig. §b). Figure||b shows that the resummed 
calculation is in better agreement with the data points than the NLO expression. We have 
found it difficult to reproduce the rapid growth of (qp) as xp — > in either approach. In 
the future, it will be interesting to see how a more precise theoretical study will be able to 
adequately explain this rapid growth of (gf,) in the region xp — > 0, assuming that the actual 
experimental data for the (qp) vs. xp distribution resemble the "derived data" discussed 
above. 

The resummation also significantly affects the pp dependence of the charged particle 
multiplicity. In Fig. [|we present the distribution (l/a tot )da/dpp. We see that resummation 
effects must be included to describe the shape of this distribution at pp < 1 GeV. Further- 
more, resummation also improves the agreement between the theory and the experiment in 
the whole range of pp. Through Eq. ( f>5"D , the improved description of the qp distribution in 
the small-gr region translates into a better agreement with the pp distribution in the whole 
range of pp. Just as in the case of the z-flow, the fixed-order calculation gives a rate that is 
too small compared to the data, which implies that higher-order corrections are important. 
Until the complete 0(a 2 s ) corrections are available, the resummation formalism, which al- 
ready accounts for the most important contributions in the region of the phase space with 
the highest rate {i.e., at small qp), serves as a better theoretical prediction in the whole 
range of pp. 

Finally, Fig. [5] shows the ^-distribution for the charged particle multiplicity 
(1/ 'cr tot )da /dxp. We see that both the resummed and fixed-order distributions are in rea- 
sonable agreement with the data and with earlier published theoretical results for the 0(a s ) 
Xir-distributions | 15fl . For the fixed-order multiplicity, we present two additional curves corre- 
sponding to different choices of the factorization scale fip in (|5]); the lower and upper dotted 
curves correspond to fip = 0.5Q and 2Q, respectively. Note that the scale dependence of 
the NLO multiplicity increases when z — > 0. Also note that the resummed multiplicity is 
significantly lower than the data in the two lowest bins of xf {(%f) = 0.075 and 0.125), but 
consistent with the NLO multiplicity within the uncertainty due to the scale dependence. 
Such behavior of the resummed multiplicity results from the dependence of the 0(a s ) co- 
efficient functions C^^ l \z, C±, C2, Hf) on the additional term lnz which was given in 
Eqs. (|72|) and (|73D and discussed at the end of Section ||. This negative logarithm dominates 
the C ou * (1) -functions at very small values of z. Similarly, the integral ([12]) of the NLO cross 
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section over the lowest bin < q\ < (q^) 2 depends on lnz through the terms 

In j^y 2 (6(1 - z)5 bj P ja (x) + P bj (z)5(l ~ ^ja) , 

as given in flI3|). Numerically, this dependence is less pronounced than in the resummed 
cross section. For z ^ 0.1, the growing scale dependence of the multiplicity in the 0(a s ) 
calculation indicates that unaccounted higher-order effects become important and are needed 
to improve the theory predictions. For example, including the 0(a 2 s ) coefficient Cj® in the 
resummed calculation will be necessary to improve the description of the charged particle 
multiplicity in the small- z region. 



V. DISCUSSION AND CONCLUSIONS 



In this paper we have continued the discussion of the effects of multiple parton radiation 
in semi- inclusive deep-inelastic scattering that was started in 0. We have demonstrated 
that these effects influence a large class of observables and can be described with the help of 
a resummation formalism M that is similar to the resummation formalism used to improve 
the theoretical description of energy-energy correlations in e + e~ hadroproduction JTIJ, and 
small-pr vector boson production at hadron-hadron colliders [jT^ . 

Although this formalism needs further improvement, in particular in the procedure for 
matching the resummed curve with the perturbative result in the transition region, it already 
improves the agreement between the theory and the data and provides interesting insights 
about qualitative features of semi-inclusive DIS. We have demonstrated that the resumma- 
tion formalism describes the pseudorapidity distributions of the transverse energy flow more 
accurately than the NLO calculation; this formalism also has good potential to improve the 
description of various distributions of particle multiplicity. 

From our analysis we conclude that the experimentally measured qx distribution for the 
energy flow broadens rapidly as x — > 0. This rapid broadening of the distribution can 
be realized if the nonperturbative Sudakov factor in the resummed energy flow increases as 
1/x. Similarly, the qx distribution for the charged particle multiplicity broadens rapidly when 
z —>■ 0, which is consistent with the nonperturbative Sudakov factor increasing as z~ lA . The 
sDIS nonperturbative Sudakov factors at small values of x and z are therefore qualitatively 
different from the known nonperturbative Sudakov factors for vector boson production and 
e + e~ hadroproduction, which do not depend on the longitudinal variables at all. The rapid 
growth of the nonperturbative Sudakov factor in sDIS might indicate that the ep collider 
HERA tests the resummation formalism in a new dynamical regime, which was not studied 
yet at colliders of other types. 

There are several theoretical aspects of the resummation formalism that can be clarified 
when more experimental data are published. Perhaps the largest uncertainty in the predic- 
tions of the resummation formalism comes from the unknown nonperturbative contributions, 
which in the 6-space formalism are included in the nonperturbative Sudakov factor S NP (b). 
In this paper, we have presented simple parametrizations of S NP (b) for the transverse energy 
flow (cf. Eq. fl35p) and charged particle multiplicity (cf. Eq. fl¥Ip). These parametrizations 
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were found by fitting the resummed energy flow and charged particle multiplicity to the data 
from Refs. ||,U and Ref. |J, respectively. Experimental measurements outside the range of 
those data will make it possible to further improve these parametrizations and, hence, the 
accuracy of the resummation formalism. 

The most straightforward way to study S NP (b) is by measuring the variation of the 
qr spectra of physical quantities due to variations of one kinematical variable, with other 
variables fixed or varying only in a small range. For the energy flow, it would be beneficial 
to obtain more data at x > 10~ 2 , where the predictions of the resummation formalism 
can be tested more reliably, without potential uncertainties due to the small-x physics. 
Another interesting question is the dependence of the nonperturbative Sudakov factor on 
the virtuality Q of the vector boson. This dependence can be tested by studying the qx 
spectra in a range of Q with sufficient experimental resolution in the current fragmentation 
region. Finally, to study effects of multiple parton radiation on semi-inclusive production 
of individual hadrons, it will be interesting to see the q? spectra for particle production 
multiplicities with the additional constraint z > 0.05 ~ 0.1, i.e., in the kinematical region 
where the parametrizations of the fragmentation functions are known reasonably well. 
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APPENDIX A: SUMMARY OF NOTATIONS 

Detailed explanations of the notations are provided in [|T| . 

We limit the discussion to the case of the semi-inclusive DIS process 

e + e + B + X (43) 

at the ep collider HERA. Here e is an electron or positron, A is a proton, B is a hadron 
observed in the final state, and X represents any other particles in the final state in the 
sense of inclusive scattering. We denote the momenta of A and B by p^ and and the 
momenta of the lepton in the initial and final states by Z M and l'^. Also, g M is the momentum 
transfer to the hadron system, q 11 = Z M — I" 1 . We assume that the initial lepton and hadron 
interact only through a single photon exchange. Therefore, q^ also has the meaning of the 
4- momentum of the exchanged virtual photon 7*. 

The kinematics of the process ( f4~3"D can be discussed conveniently with the help of the 
following five Lorentz scalars: 

• the square of the center-of-mass energy y/S e A of the initial hadron and lepton, 

SeA ={PA + If = 2 PA ■ k (44) 
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• the momentum transfer from the leptonic to the hadronic part of the process, 

Q 2 = - q 2 = 11 ■ (45) 

• the Bjorken variable, 

x = -^—; (46) 
2pa • q 

• a variable z describing energy transfer from the composite photon-proton intermediate 
state to the observed final hadron B, 



_ Pb-Pa _ 2xPb ■ Pa , (A7 , 
q-PA Q 2 



• a scale qx that serves as a control parameter in the resummation formalism, 

q 2 T = -fa*, (48) 

where 

a a n q-PB n q-PA fAn . 

<K = q^ - Pa Pb • (49) 

Pa -Pb Pa- Pb 

Let a denote the parton in A that participates in the hard scattering, with momentum 
Pa = CaP% Let b denote the parton of which B is a fragment, with momentum = p B /£b- 
We use the symbol " " to denote quantities at the parton level. Correspondingly, we 
introduce the parton-level Lorentz scalars x, z that are analogous to the hadron-level scalars 

x, z: 

x = = -, (50) 

^ Pb-Pa Z 

z= = — . 51 

q-Pa €b 

The paper concentrates on the discussion of experimental data presented in the hadronic 
center-of-mass (hCM) coordinate system, which is defined to be the center-of-mass frame 
of the proton A and the virtual photon 7*, with the z axis pointing in the direction of the 
photon's momentum. The momenta of A and 7* in the hCM frame are 

A = (^,0,0,-^1, ( M ) 



2W 2W 
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where 

W = Qj- - 1. (54) 
V x 

The momentum of the outgoing hadron B in this frame lies in the xz plane, and its component 
perpendicular to the z axis is 

p T = zq T . (55) 
The energy and longitudinal component of the momentum of B are 

Eb = z-^-, (56) 
W " ~ It 

Pr = z — • (57) 

The pseudorapidity of B in this frame is simply related to qx, by 

qT = We~ Vcm . (58) 

The region of small q?, where our resummation calculation is relevant, corresponds to the 
region of large pseudorapidity in the hCM frame. 

We denote by (p the angle between the xz plane (the plane of the momenta pa and 
Pb) and the plane of the lepton momenta i and £'. In terms of the angle p and the boost 
parameter ip = cosh" 1 (2xS £ a/Q 2 — 1), the momenta of the leptons in the hCM frame are 

F = \-^-((W 2 + Q 2 )coshiP + l^ 2 -Q 2 Y-sinh^cos^,--sinh^siny?, 



\4W 
1 

w 



(W 2 + Q 2 ) + (W 2 -Q 2 )coship)y, (59) 



= i JLf(^ 2 + Q 2 )cosh^- W /2 + Q 2 ),^sinh^cos^,-^sinh^sin^, 

1 ' w 2 -g 2 + (i^ 2 -g 2 )cosh^U. (60) 



AW 

It is convenient to decompose the sDIS cross section into a sum over a complete set of 
functions A p (tp, (p) of the azimuthal angle <p and the parameter ip [|7|: 

1 jE^ = ^v B ^M\&)A^M- (6i) 

If the scattering of leptons on hadrons is caused by a parity-conserving interaction, then 
the spin-averaged sDIS cross section receives contributions only from four angular functions 
A p (if),ip), which are 

A 1 = 1 + cosh 2 ip, 

A 2 = -2, 

A 3 = — cospsmh2il), 

At = cos2<£> sinh 2 ip. (62) 
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Since the angular functions A 3 and A4 are proportional to cos if and cos 2y?, respectively, they 
do not contribute to the cross section that is integrated over the azimuthal angle (p. Hence, 
the (^-integrated cross section da / (dxdzdQ 2 dq^) depends only on the structure functions Ai 
and A 2 . 



APPENDIX B. PERTURB ATIVE COEFFICIENTS IN THE RESUMMED CROSS 

SECTION 

The Wilson coefficient functions and C£"*, and the perturbative part S p (b, Q) of the 
Sudakov factor in (^0) can be expanded in powers of a s WI]], as 

2n2 



CiQ 



V 2 



JCf/b 2 /J, k=l \ 7T J L 

+ Bf^a^pj.d,^)}; (63) 

cz($,c u c 3 ,b,n F ) = EC w (^^)f^^V; (64) 

fc=0 V t / 

C£*{z,C u C a ,b,» F ) = ±C^ k \z,^b)(^^) k . (65) 



fc=0 



Up to the order 0(a s ), the coefficients in Eqs. (|63])- (|65|) are listed below. 

• For the Sudakov factor, 

A {1) = C F , (66) 

Here C F = 4/3 is a color factor; fj, F is a factorization scale; 60 = 2e _7E ; C\ and C2 are 
constants that are conventionally chosen such that C1/C2 = &o- 

• For the coefficient functions, 

Cf°\i^ F b) = = M(i - 0, (68) 

C M0) = C ou m = Q) (69) 

Clf 1] (x,fx F b) = - x) - P qq {x)\n(^- 

- w( i-fi)(5 + m«(^)), (70) 

C^\i it iFb) = \x{l-x)-P qg {x)\n(^-), (71) 
C^W) = ^(l-2)-P„(2)ln(^) 
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,16 V fe C7 2 

C^ {1) = C°f (1) = 0. (74) 

Here P gg and P ?ff are one-loop splitting functions 0. The subscripts a, b correspond 
to quarks, antiquarks and gluons; the subscripts i,j correspond to the quarks and 
antiquarks only. 
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FIG. 1. Comparison of the resummed z-flow in the current region of the hCM frame with the 
data in the low-Q 2 bins from Refs. Q (filled circles) and || (empty circles). The parametrization 
of the nonperturbative Sudakov factor in (^) is used. For the bin with (Q 2 ) = 33.2 GeV 2 and 
(x) = 0.0047, the 0(a s ) contribution for \xf = Q is shown as a dashed curve. 
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FIG. 2. Comparison of the resummed z-flow in the current region of the hCM frame with the 
data in the high-Q 2 bins from Ref. j| . The parametrization of the nonperturbative Sudakov factor 
in (|35|) is used. For the bin with (Q 2 ) = 617 GeV 2 and (x) = 0.026, the 0(a s ) contribution for 
fJ-F = Q is shown as a dashed curve. 
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FIG. 3. The distributions (a) {pj) vs. xp and (b) (g|,) vs. xp for the charged particle multiplic- 
ity at (W) = 120 GeV, (Q 2 ) = 28 GeV 2 . The experimental points for the distribution (pj,) vs. xf 
are from Fig. 3c of Ref. ||. The "experimental" points for the distribution {qj) vs. xf are derived 
using Eq. (f42|). The solid and dashed curves correspond to the resummed and the next-to-leading 
order (/ip = Q) multiplicity, respectively. 
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FIG. 4. The dependence of the charged particle multiplicity on the transverse momentum pt 
of the observed particles in the hCM frame. The data points are from ||. The solid and dashed 
curves correspond to the resummed and NLO multiplicities, respectively. 
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FIG. 5. The dependence of the charged particle multiplicity on the Feynman variable xf in 
the hCM frame. The solid curve corresponds to the resummed multiplicity. The dashed, lower 
dotted and upper dotted curves correspond to the NLO multiplicity calculated for \xf = Q, 0.5Q 
and 2Q, respectively. 
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